Multistability at arbitrary low optical intensities in a metallo-dielectric layered 

structure 
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We show that a nonlinear metallo-dielectric layered slab of subwavelength thickness and very 
small average dielectric permittivity displays optical multistable behavior at arbitrary low optical 
intensities. This is due to the fact that, in the presence of the small linear permittivity, one of the 
multiple electromagnetic slab states exists no matter how small is the transmitted optical intensity. 
We prove that multiple states at ultra-low optical intensities can be reached only by simultaneously 
operating on the incident optical intensity and incidence angle. By performing full wave simulations, 
we prove that the predicted phenomenology is feasible and very robust. 

PACS numbers: 78.67.Pt, 42.65.Tg 



Hysteresis is probably one of the most intriguing fea- 
ture a nonlinear system can exhibit both theoretically 
and for its main application, i.e. designing of memory 
devices. In optics, bistability and related hysteresis have 
attracted a large research interest in the last decades 0, 0] 
since they are the basic ingredients for devising optical 
memories, logic gates and optical computing devices 
However, the actual exploitation of the hysteresis behav- 
ior has been hampered by the fact that the electromag- 
netic bistable phenomenology generally occurs at high 
optical intensities. The research effort aimed at reduc- 
ing the bistability threshold is considerable and almost 
all considered strategies rely on the enhancement of non- 
linear effects [3-0] • Metallo-dielectric multilayer struc- 
tures have been shown to host considerable nonlinearity 
enhancement p, 3| and this has lead to the prediction 
of bistability jloL characterized by a low intensity 
threshold Optical bistability at relatively low opti- 
cal intensities has also been predicted in the presence of 
surface plasmon polaritons |13| where the low bistability 




threshold is due to the strong plasmon field enhancement 
which in turn produces an efficient enhancement of the 
nonlinear effects. Following a very different route, it has 
recently been proposed that an extremely marked non- 
linear behavior can be observed, instead of by enhanc- 
ing the nonlinear response, by substantially reducing the 
medium linear dielectric permittivity [l4| and this strat- 
egy has allowed to predict low-intensity transmissivity 
directional hysteresis exhibited by nonlinear metamate- 
rial slabs with very small linear permittivity [l5j . 

In this Letter, we show that a nonlinear metallo- 
dielectric layered slab of subwavelength thickness and 
very small average dielectric permittivity is character- 
ized by a multistable behavior which is accessible at ar- 
bitrary low optical intensities. More precisely, we prove 
that the slab transmissivity, in addition to the angular 
multistable behavior occurring at low (but not arbitrary 
small) optical intensities (as discussed in Ref. 15]), is a 
multi-valued function of the incident optical intensity at 



FIG. 1: (Color online) Geometry of the layered metallo- 
dielectric slab and of the TM incident (i), reflected (r) and 
transmitted (t) plane waves. 



each fixed incidence angle. The nonlinear electromag- 
netic matching at the output free-space slab interface 
combined with the small linear permittivity produces the 
multiplicity of electromagnetic states, one of which has 
a large electric field component normal to the slab even 
if the free-space transmitted intensity approaches zero. 
As a consequence, the lowest intensity at which trans- 
missivity is multi- valued turns out to vanish for normal 
incidence and this implies that, at a very small incidence 
angle, multistability occurs at ultra-low optical intensi- 
ties. The overall slab transmissivity turns out to be a 
rather articulated multi- valued function of both incident 
optical intensity and incidence angle and we show that 
ultra-low intensity multistable states can be reached only 
by simultaneously varying the incident optical intensity 
and incidence angle along suitable paths. We perform full 
wave simulations to investigate both the nontrivial meta- 
material homogenization in the strong nonlinear regime 
and the role played by plasmonic resonances. We show 
that the induced plasmon polaritons excited at the metal- 
dielectric interfaces affect the numerical value of the slab 
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FIG. 2: (Color online) Nonlinear slab transmissivity T (solid 
line) as a function of the normalized input field intensity Ii„ 
at the fixed incident angles 6 = 0.028 rad (panel (a)) and 
9 = 0.170 rad (panel (b)). (c) Surface \x\ 1/2 A z and values of 
|x| 1,/2j 4z (solid lines) corresponding to the transmissivities of 
panels (a) and (b). In the three panels, capital letters label 
some reference states. 



transmissivity solely in proximity of the hysteresis jumps 
without altering the overall hysteresis behavior which is 
therefore proven to be a feasible and robust phenomenol- 
ogy. 

Consider the metallo-dielectric layered slab reported in 
Fig.l, embedded in vacuum, of thickness L = 258 nm and 
consisting of alternating, along the y-axis, two isotropic 
and non-magnetic nonlinear Kerr layers of thicknesses 
di and c?2- We choose to operate with a monochro- 
matic radiation of wavelength A = 810 nm and, if 
the layers thicknesses are such that d\ + d,2 <C A, the 
slab behaves as an homogeneous medium which, in the 
presence of a transverse magnetic field (TM), is char- 
acterized by the standard Kerr constitutive relations 
D = e {eE + x [(E • E*) E + \ (E • E) E*] } [14]. The 
effective dielectric permittivity e and nonlinear suscepti- 
bility x are the average values of their underlying con- 
stituents and we choose here e = —0.054 and \ — 
8.72 x 10~ 17 m 2 /V 2 (see below for the actual feasibility 
of these values). Note that the negative and very small 
value of the effective dielectric permittivity is achieved 
here by averaging the negative and positive permittiv- 
ities characterizing metal and dielectric layers, respec- 
tively. As reported in Fig.l, a TM plane wave of am- 
plitude Ei is made to impinge onto the slab interface 
at z = with incidence angle 9 thus producing a TM 
reflected plane wave of amplitude E r and a transmit- 
ted TM plane wave of amplitude E t . We have numeri- 
cally evaluated the slab transmissivity T = \E t /Ei\ 2 for 



a number of incident optical intensities and incidence an- 
gles and in Fig. 2 we report T as a function of the input 
intensity 7j„ = (l/2)y / eo/ u \Ei\ 2 for the incidence an- 
gles 6 = 0.028 rad (panel (a)) and 6 = 0.170 rad (panel 
(b)). The reported slab transmissivities are multi- valued 
functions of the incident optical intensity. The multiplic- 
ity (at last three-fold) of different slab electromagnetic 
states, corresponding to the same transmissivity, phys- 
ically arises (as discussed in Ref.[ll|) from the impact 
of the nonlinearity on the electromagnetic matching con- 
dition at the slab interfaces, an effect made possible by 
the very small value of |e|. Specifically, the continuity, 
at z = L, of the displacement field component normal 
to the output slab interface, D Z (L~) = D Z (L + ). yields 
[e + lx{\E t \ 2 cos 2 B + Afj] A z + \Et\sm9 = 0, which is 
a cubic equation that, for a given E t and yields at 
last three different values of the real amplitude A z = 
(Et/\Et\) exp(— icos9)E z (L) characterizing the longitu- 
dinal electric field component at z = L. In Fig. 2(c) 
we report the surface comprising all the possible values 
of \x\ x ^ 2 A z together with values of \x\ x ^ 2 A z (solid lines) 
corresponding to the transmissivities of Fig. 2(a) and (b), 
the capital letters labelling some reference states. Note 
that the incident intensities at which multistability oc- 
curs (of the order of 0.15 MW/cm 2 and 4 MW/cm 2 for 
the transmissivities of Fig. 2(a) and 2(b), respectively) are 
smaller than those normally required for observing stan- 
dard optical bistability. However it is evident that the 
smaller |e| the smaller the intensity threshold required 
for observing the considered multistability so that very 
smaller intensity thresholds (of the order of W/cm 2 ) are 
very likely to be attained. 

Multistable transmissivities of Fig. 2(a) and 2(b) show 
two main mutual differences which deserve to be dis- 
cussed. First, multi- valued transmissivity of Fig. 2 (a) has, 
for Ii n > In, three different branches that does not break 
whereas multi-valued transmissivity of Fig. 2(b) has, for 
Iin > /_e, an upper branch that does not break and two 
lower branches characterized by the two breaking states 
B and H. The origin of this difference is easily under- 
stood by considering the surface of Fig. 2(c) from which 
it is evident that the breaking of a transmissivity branch 
occurs for states close to critical points at which surface 
folding takes place (states B and H of Fig. 2(c)). Second, 
even though in both cases there exists a lower intensity 
above which the transmissivity is multi-valued (see the 
state N and E of Fig. 2(a) and 2(b), respectively), trans- 
missivities of Fig. 2(a) and 2(b) are multistable for I in > 
I N = 0.04 MW/cm 2 and I in > I E = 1.92 MW/cm 2 , re- 
spectively so that Ijv Ie or, in other words the smaller 
the incidence angle the smaller the incident optical in- 
tensity at which multistable states are accessible. In or- 
der to closer investigate this remarkable point we have 
plotted in Fig. 3 the slab transmissivity as a function of 
both the incident optical intensity and incidence angle 
thus providing a complete description of the slab mul- 
tistability. The transmissivity of Fig. 3 is a multi-sheet 
surface along which we have also plotted, for clarity pur- 
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FIG. 3: (Color online) Slab transmissivity T as a function of 
both the incident optical intensity Ii„ and incidence angle 6. 
Dashed lines represent the slab transmissivity of Fig.2(a) and 
2(b) and the capital letters label the same reference states as 
in Fig.2. In the inset the plane of independent parameters 
Tin and 9 is reported together with the region at which slab 
multistability occurs (shaded region). 



poses, the transmissivities of Fig. 2(a) and 2(b) (dashed 
lines) together with the very same reference states la- 
belled with capital letters. The intersection between the 
upper and intermediate sheets yields the lower intensity 
at which multistability occurs and it turns out that such 
intensity is a monothonically increasing function of inci- 
dence angle 9 and, most importantly, that this intensity 
vanishes if the angle vanishes. To underline this result, 
in the inset of Fig. 3 we have reported the region of the 
independent external parameter plane , 9) where mul- 
tistability occurs (shaded region) and, since this region 
contains the excitation state 7j n = 9 = 0, it is evident 
that multistable electromagnetic slab states are available 
at arbitrary small incident optical intensities (and very 
small incidence angles). In order to physically grasp the 
origin of such an unusual situation, consider the states 
providing multistability and belonging to the lower sheet 
of the surface reported in Fig. 2(c). For |i£ t | = we have 
A z = and A z = ±- v /-2e/(3x) = ±0.19 /y/x (see the 
above mentioned matching condition at z = L) so that 
the surface exists even at arbitrary small values of \E t \ 
(see Fig. 2(c)) and multistable states exist even in the lim- 
iting case where there is no transmitted waves. Since the 
input intensity is smaller than or equal to the output one, 
it is evident that the described mechanism provides mul- 
tistable states even at ultra-low input optical intensities. 

The transmissivity reported in Fig. 2(b) clearly shows 
hysteresis behavior since the state B is a breaking point 
of its transmissivity branch. Suppose to switch on the 
incident plane wave with 9 = 0.17 rad and Ii n — I a — 
0m5MW/cm 2 so that the state A of Fig.2(b), 2(c) and 
Fig. 3 is excited since before the illumination the slab was 
not polarized. By increasing the input intensity and hold- 
ing the incidence angle fixed, the transmissivity follows 



the curve joining the point A and B of Fig. 2(b) while the 
electromagnetic state continuously varies as in Fig. 2(c). 
From the state B, if the intensity is further increased, 
it is evident that the electromagnetic state undergoes a 
sudden jump to the state C (see Fig. 2(c)) on the lower 
surface sheet since there is no allowed state continuously 
joined to B. As a consequence, the transmissivity under- 
goes a sudden jump to a higher value (see Fig. 2(b)). If 
now, starting from the state C, the input intensity is de- 
creased, the electromagnetic state varies along the curve 
from C to D of Fig. 2(c) whereas the transmissivity in 
Fig. 2(b) assumes the values along the curve from C to 
D which arc different from those attained along the for- 
ward path (hysteresis). At the state D, if the intensity 
is further decreased, the state undergoes a jump to the 
state E of Fig. 2(c), belonging to the central surface sheet, 
since the states from D to F of Fig. 2(c) have, according 
to Fig. 2(b), higher input intensities. 

The transmissivity reported if Fig. 2(a), although 
multi-valued, does not exhibit hysteresis since its 
branches do not present breaking points. However the 
low intensity multistable states can always be reached 
by selecting appropriate paths of the excitation plane 
(Iin, 9). Suppose, for example, the multistability respon- 
sible state L of Fig. 2(a), 2(c) and Fig. 3 (of intensity 
II = 0.08 MW/cm 2 ) is required to be reached. A possi- 
ble operational way consists in exciting the state A and 
reaching the state C as discussed in the above hystere- 
sis example. Now if from the state C the intensity is 
hold fixed and the angle is decreased, from Fig. 3 it is 
evident that the state I is reached since both the states 
C and I belongs to the same surface sheet of Fig. 2(c). 
From the state C one can now fix the angle and decrease 
the incident optical intensity thus eventually reaching the 
prescribed target state L. From L one can decrease the 
intensity thus reaching the state M from which, a fur- 
ther intensity decrement, produces (in analogy to the 
above described hysteresis loop) a sudden jump of the 
trasmittivity to N of Fig. 2(a), the corresponding state 
N of Fig. 2(c) belonging to the central surface sheet. The 
very same procedure can even be used to excite the above 
discussed multistable states at much smaller incident op- 
tical intensities (and incidence angle). 

In order to discuss the feasibility of the above predicted 
slab multistability we have performed 3D full- wave finite- 
element simulations 16] for evaluating the transmissivity 
of the slab in the presence of a TM radiation of free- 
space wavelength A = 810 nm. We have chosen d% = 2 
nm, d 2 = 5.25 nm, L = 258 nm, ei = -28.81 + lOi, 
e 2 = 10.9 - 3.8i, xi = 3.16 x 10~ 16 m 2 /V 2 and X 2 = 0. 
The parameters of medium 1 coincide with those of sil- 
ver [17] . characterized by a very large nonlinear suscep- 
tibility [l8j . with the imaginary part of the permittiv- 
ity corrected by the layer size effect (since c?i = 2 nm) 
[13 ] whereas medium 2 is a linear dielectric with gain (to 
compensate the metal losses). The effective permittiv- 
ity and nonlinear susceptibility of the considered sample 
are e = (did + <2 2 e 2 )/(di + d %) = -0-054 + 0.007? and 



4 



!-5(a) 




FIG. 4: (Color online) (a) Comparison between the slab trans- 
missivities evaluated through full-wave simulations (dotted 
lines) and those of Fig. 2(b) (solid lines), (b) Semi-log plot 
of the maximum values, within the layered medium bulk, of 
the normalized squared field components as a function of the 
normalized input field amplitude obtained for the full-wave 
evaluation of the slab transmissivity reported in panel (a). 



X = (diXi + <feX2)/(di + da) = 8.72 x lfr 17 m 2 /V 2 . In 
Fig. 4(a) we report the comparison between the transmis- 
sivity evaluated through full-wave simulations for 8 = 
0.17 rad (dotted line) and the transmissivity of Fig. 2(b). 
We note that good agreement exists between the results 
of the two kind of simulations and, most importantly, 
that finite-element simulations still predicts the above 
discussed multistability thus proving its robustness. The 
origin of the discrepancies between the transmissivities 
lies in the fact that a surface plasmon resonance occurs 



when the TM wave impinges on the slab as reported in 
Fig.l. Therefore an electric field y-component (E y ) of 
plasmonic origin arises mainly at the edges of the lay- 
ers inside the medium and it is characterized by a sub- 
wavelength varying profile and an evanescent field tail 
in vacuum which does not contribute to the power flow. 
Evidently, if E y is much smaller than E x and E Zl the 
homogenization theory correctly describes the slab non- 
linear behavior. In Fig. 4(b) we report the maximum 
values (within the medium) attained by the three field 
components for 8 = 0.17 rad as a function of the in- 
cident optical intensity corresponding to the evaluated 
full-wave transmissivities reported in Fig. 4(a). It is ev- 
ident that, outside the shaded region, E y (dot-dashed 
line) is much smaller than both E x and E z (solid and 
dashed lines respectively) and in the corresponding re- 
gions of Fig. 4(a) the agreement with the homogenization 
approach is very satisfactory. Within the shaded region 
of Fig. 4(b), E y is not negligible and this partially breaks 
the validity of the homogenization approach thus leading 
to the discrepancies of the transmissivities in the range 
1 MW/cm 2 < I in < 4 MW/cm 2 as reported in Fig.4(a). 

In conclusion we have shown that a nonlinear metallo- 
dielectric layered slab of subwavelength thickness with 
very small average permittivity is characterized by a very 
low bistability threshold and, most importantly, that it 
provides multistable states at arbitrary low optical in- 
tensities. Such ultra-low intensity multistable states are 
always accessible through a suitable excitation procedure 
involving both the incident intensity and the incidence 
angle. Since the main difficulty of usefully exploiting 
standard optical bistability is related to the required high 
values of the optical intensity, we believe that our findings 
can have a large impact on the devising of a novel class 
of nanophotonic devices with complex memory function- 
alities and operating at very low optical intensities (as 
compared to those generally required in standard nonlin- 
ear optics). 
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